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AMENDMENTS TO THE CLAIMS 

1. (CURRENTLY AMENDED) A method of interactively simulating 
contact between a deformable first object (201; — 2^-; — 4r0-; — &9-j — 
and a second object (202 ; — ^ — — ^9-? — £0-? — 91) using a simulated 
model with a predetermined sampling time step, the method being 
characterized in that: 

(a) the parameters describing the physical characteristics 
of each of the objects, such as the geometry and the mechanics 
of the materials of each of the objects, are calculated 
beforehand and stored in a memory, 

(b) at the beginning of each sampling time step of the 
simulated model, a real-time analysis of the inherent behavior 
of each object is carried out at the level of each object in 
order to predict the positions, speeds and accelerations of that 
object in application of a free movement that does not take 
account of any subsequent contacts, 

(c) in each sampling time step of the simulated model, 
pairs of objects that are detected as intersecting are analyzed 
in real time at the level of an overall scene including the 
objects liable to come into contact, and a list of groups of 
collisions is established that contains a string of objects in 
collision and a description of the collisions, 

(d) in each sampling time step of the simulated model, 
parameters representing the physical characteristics of the 
objects and the description of the collisions are repatriated in 
real time for each group of collisions to determine, for each 
instance, the solution to the Signorini problem that governs 
contact between two objects in the case of pure relative 
sliding, 

-4- 

WEINGARTEN, SCHURGIN, 
GAGNEBIN & LEBOVICI LLP 
TEL. (617) 542-2290 
PAX. (617) 451-0313 



Attorney Docket No. BDL-506XX 
Filed: Herewith 
TC Art Unit: 
Confirmation No. : 

(e) at the end of each sampling time step of the simulated 
model, a real-time display of the inherent behavior of the 
object following the collision is effected at the level of each 
object, and 

(f) all real-time processing is effected with a 
computation time step shorter than the sampling time step of the 
simulated model so as to define an interactive simulation in 
which the user can intervene directly during simulation. 

2. (ORIGINAL) A method according to claim 1, characterized in 
that during the step a) of calculating beforehand parameters 
describing the physical characteristics of each of the objects, 
a finite element type description of deformations is used for 
the parameters describing the mechanics of the materials, with 
matrices being filled and inverted, systems of equations being 
solved, and data being stored in memory. 

3. (CURRENTLY AMENDED) A method according to claim 1 or claim 
2, characterized in that each object is described in a rest 
configuration as a set of triangles reproducing its surface and 
a set of tetrahedra describing the interior of the object. 

4. (ORIGINAL) A method according to claim 3, characterized in 
that each triangle is described by three points placed in an 
order such that normals are calculated that are invariably 
directed towards the exterior of the object. 

5. (CURRENTLY AMENDED) A method according to claim 3 or claim 
4, characterized in that the deformations of the objects are 
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interpolated by the finite element method using a linear 
tetrahedral mesh. 

6. (CURRENTLY AMENDED) A method according to any — ene — e# 
claims- 1— fee — 5-, characterized in that in each computation time 
step the explicit forces applied to an object, which are already 
known at the start of the computation step, are integrated 
during the step b) at object level to define the movements of 
the object that they generate, whereas the values of the 
implicit contact forces, which depend on the movement of the 
objects in the computation time step, are determined during the 
step d) of seeking the solution to the Signorini problem at the 
level of an overall scene. 

7. (CURRENTLY AMENDED) A method according to any — eae — e# 
claims 1 — fee — &, characterized in that during the step c) of 
analysis at the level of an overall scene, the existing 
intersections between the objects of the scene are detected 
geometrically in order to extract from pairs of elements of 
intersecting objects a length and a direction of 
interpenetration between the two elements of a pair of elements 
of objects. 

8. (ORIGINAL) A method according to claim 7, characterized in 
that during the step c) of analysis at the level of an overall 
scene, to extract from pairs of elements of intersecting objects 
a length and a direction of interpenetration between the two 
elements of a pair of elements of objects, an intermediate 
movement of the objects between the preceding computation step 
and the current computation step is also taken into account in 
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order to compute a preferential direction of interference 
between the objects. 

9. (CURRENTLY AMENDED) A method according to afty — one — e# 
claimo — 3 — fee — — and according to claim 7 claim 24 , characterized 
in that during the step d) of seeking the solution to the 
Signorini problem, the extreme points of application of the 
contact force between the two objects in collision are 
reconstructed if those extreme application points have not been 
determined during the preceding step. 

10. (CURRENTLY AMENDED) A method according to claim 9, 
characterized in that during the step d) of seeking the solution 
to the Signorini problem, in the case of a segment -segment 
intersection of two triangular objects — — &Of f the two points 

( 4 1, 5^ selected to constitute the extreme points of 

application of the contact force between the two objects (40 , 
5-0-) — in collision are situated at the intersection of each of the 
two segments (PiP 2 , Q1Q2) with the plane formed by the face of the 
triangle in the intersection. 

11. (CURRENTLY AMENDED) A method according to claim 9, 
characterized in that during the step d) of seeking the solution 
to the Signorini algorithm, in the case of a point- face 
intersection of two triangular objects — f^Oi — Wf , a first point 
(71) — selected to constitute a extreme point of application of 
the contact force between the two objects (60, — 7^) — in collision 
is the point of the intersection whereas the second extreme 
point of application of the contact force between the two 
objects in collision is the projection ( 6 1) — of the first extreme 
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point 4^ onto the face of the triangle 46^) in the 

intersection. 

12. (CURRENTLY AMENDED) A method according to any — ene — e# 
claime 9 — fee — ii, characterized in that barycentric coordinates 
are used to distribute the displacements and the forces of the 
points of application of the contact force between the extreme 
points of application of the contact force by effecting a linear 
interpolation for a finite element modeling process. 

13. (CURRENTLY AMENDED) A method according to claim 25 a 10 and 
characterized in that the distance 8 of interpenetration 

between the two extreme points 443^ — — of application of the 
contact force in the case of a segment -segment contact between a 
first segment (Qi Q 2 ) and a second segment (Pi P 2 ) of a second 
triangle is calculated from the following equation: 



5 = [a s b { cj 



[a 1 - a] 



W 2J 



- [P 1-0] 



v 2 



(1) 



in which: 

a and 1 - a are the barycentric coordinates on the first 
segment (Qi Q 2) 

P and 1 - P are the barycentric coordinates on the second 
segment (Pi P 2 ) , 

ai bi ci are the coordinates of the interpenetration 
direction n it 

Wi and W 2 are the coordinates of the first segment Qi Q 2/ and 
Vi and V 2 are the coordinates of the second segment Pi P 2 . 
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14. (CURRENTLY AMENDED) A method according to claim 26 o 11 and 
characterized in that the distance 6 of interpenetration 
between the two extreme points -f&^b — ^ — of application of the 
contact force in the case of a point -plane contact between a 
point (71) of a second triangle and a plane (Pi P 2 P3) of a first 
triangle is calculated from the following equation: 











8 = [a, b, cj 


[a 0 y] 


w 2 


- v, 






w 3 





in which: 



a, P and y are the barycentric coordinates on the first 
triangle, 

ai bi ci are the coordinates of the interpenetration 
direction n ix 

Wi, W 2 , W 3 are the coordinates of the first triangle 
(Pi P 2 P 3 ), 

Vi represents the coordinates of the point of contact 
consisting of a vertex (Qi) of the second triangle (Qi Q 2 Q3) • 

15. (CURRENTLY AMENDED) A method according to any — ene — a# 
claims 1 — fee — ±4r t characterized in that when the points of 
application of the contact forces between two objects in 
collision have been determined, during the step d) the 
mechanical characteristics of the objects are transferred into 
the defined contact space in which the whole of a group of m 
contacts with n objects is processed, where m and n are 
integers . 

16. (ORIGINAL) A method according to claim 15, characterized 
in that during the step d) the mass and inertia of an object are 
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considered lumped together at its centre of mass and an 
instantaneous relationship between the contact forces f c in the 
contact direction, the accelerations 5" c caused by the 
constraints in the same direction, and the free accelerations 
5" free in the same direction known during the step c) at the level 
of an overall scene is established from the equation: 

(5" c = JcM-Vfc+^'ftec (3) 

in which: 

J c is an m*6n Jacobean matrix that transfers the 
instantaneous linear and angular movement into the contact 
space, 

J C T is the transposed matrix of J c/ 

M is a block diagonal matrix corresponding to the mass and 
inertia of the n objects of the group of contacts. 

17. (CURRENTLY AMENDED) A method according to any — eee — a# 
claimo — 9 — fee — 14 — — cither — a€ — claima — 3r5 — and — 3r& claim 27 , 
characterized in that during the step d) , to transport the local 
mechanical characteristics, a relationship is established 
between: 

the displacement difference (Uk 1 ) of the points of the 
deformable mesh representing the object i at the time k, between 
the free deformation (uVfree) and the constrained deformation 
(uVc) , in other words XJ k x = U\ iC - UVfree 

the free and constrained relative positions 5f ree and 8 C of 
the objects in the contact space: 

5 = Si=i n N c j U k j + 6 free (4) 
where Nc 1 is a matrix for passing from the displacement space of 
the mesh to the displacement space of the contacts, and 
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a relationship is established between the forces in the 
contact space f c and the forces in the deformation forces space 
F k : 

F k = (N c j ) T / c (5) 

18. (ORIGINAL) A method according to claim 1, characterized in 
that during the step d) an instantaneous linear relationship 
characterizing contact deformations or displacements 5 C from the 
contact forces f c and the free displacements 5" fre e caused by free 
movements integrating only the forces known explicitly at the 
beginning of the computation time step is established from the 
following equation: 

S c = [z M n N; A (U k _.) (Ne 1 ) 1 ] / e + bfree ( 6 ) 
in which: 

is a matrix for passing from the displacement space of 
the mesh to the displacement space of the contacts, 
(N X C ) T is the transposed matrix of N x c , 

A is a matrix defining the deformation of the object at the 
local level, such that if U k represents the vector of the 
displacement in the local frame of reference of the object at 
the current time and U k -i represents the displacement vector in 
the local frame of reference of the object in the preceding 
calculation step, the instantaneous values whereof are known at 
the beginning of the current computation step, then: 
U K = A (U k -i) F k + b (Uk-i) (7) 
where : 

F k is a vector representing the external forces applied to 
the object expressed in the local frame of reference, and 
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b is a vector that has a value in the displacement space 
and depends on the object deformation model. 

19. (ORIGINAL) A method according to claim 17, characterized 
in that during the step d) an instantaneous relationship 
characterizing the contact deformations or displacements 5 C from 
the contact forces f c and the free displacements 5" free caused by- 
free movements integrating only the forces known explicitly at 
the beginning of the computation time step is established from 
the following equation: 

S c = [e dt 2 J c M" 1 J C T +E i=1 n N c j A (U k .,) (N c l ) T ] / c + Sfree (*) 
in which: 

J c is an m*6n Jacobean matrix that transfers the 
instantaneous linear and angular movement into the contact 
space, 

J C T is the transposed matrix of J c , 

M is a block diagonal matrix corresponding to the mass and 
inertia of the n objects of the group of contacts, 

0 is a constant depending on the time integration method, 

N L C is a matrix for passing from the displacement space of 
the mesh to the displacement space of the contacts, 

(N^) 7 is the transposed matrix of N^, 
A is a matrix defining the deformation of the object at the 
local level such that if U k represents the vector of the 
displacement in the local frame of reference of the object at 
the current time and U k -i represents the displacement vector in 
the local frame of reference of the object in the preceding 
calculation step the instantaneous values whereof are known at 
the beginning of the current computation step, then: 
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U K = A (Uk-i) F k + b (U k -i) (7) 
where : 

Fk is a vector representing the external forces applied to 
the object expressed in the local frame of reference, and 

b is a vector that has a value in the displacement space 
and depends on the object deformation model. 

20. (CURRENTLY AMENDED ) A method according to any — one — e# 
claime 1 — fee — 3r£, characterized in that it further comprises a 
step of coupling with a haptic interface module to produce 
haptic sensation feedback to a mechanical system by means of 
which an operator manipulates the objects in a virtual scene. 

21. (CURRENTLY AMENDED) A system for the interactively 
simulating contact between a deformable first object (201; — 2r&r 
10; 60; 8 0) and a second object (202; 30; 50; — — 9G-; — 91) using 
a simulated model with a predetermined sampling time step, the 
system being characterized in that it comprises: 

(a) a module (100) — for calculating beforehand parameters 
describing the physical characteristics of each of the objects, 
such as the geometry and the mechanics of the materials of each 
of the objects, 

(b) a memory (102) for storing parameters previously 

calculated in the computation module — (100) , 

(c) a coupling module (101) — to a user interface (104) 
comprising a mechanical system held by a user to exert virtual 
forces on said objects in a scene of the simulated model, 

(d) a display screen (107) — for displaying said objects 
represented in the form of meshes, 
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(e) a central processor unit (100) — associated with input 
means — (103) , comprising at least: 

el) an object analysis module for analyzing in real 
time at the level of each object the inherent behavior of 
the object in order to predict the positions, speeds and 
accelerations of that object in application of a free 
movement that does not take account of any subsequent 
contacts, 

e2) an analysis module for an overall scene including 
the objects liable to come into contact, for analyzing in 
real time pairs of objects that are detected to be 
interacting and to establish a list of groups of collisions 
that contains a string of objects in collision and a 
description of the collisions, 

e3) a module for the real time repatriation, for each 
group of collisions, of the parameters representing the 
physical characteristics of the objects and the description 
of the collisions, for determining, for each instance, the 
solution to the Signorini problem that governs contact 
between two objects in the case of pure relative sliding, 

e4) a module for processing each object for real time 
display at the level of each object of the inherent 
behavior of that object following a collision, and 

e5) means for determining a computation step shorter 
than the sampling time step of the simulated model so as to 
define an interactive simulation. 

22. (CURRENTLY AMENDED) A system according to claim 21, 
characterized in that it comprises means for producing haptic 
sensation feedback to the user interface — (104) . 
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23. (CURRENTLY AMENDED) A system according to claim 21 — 
claim 22 , characterized in that the computation step corresponds 
to a frequency greater than or equal to approximately 500 Hz. 

24. (NEW) A method according to claim 3, characterized in that 
during the step c) of analysis at the level of an overall scene, 
the existing intersections between the objects of the scene are 
detected geometrically in order to extract from pairs of 
elements of intersecting objects a length and a direction of 
interpenetration between the two elements of a pair of elements 
of objects. 

25. (NEW) A method according to claim 10, characterized in 
that barycentric coordinates are used to distribute the 
displacements and the forces of the points of application of the 
contact force between the extreme points of application of the 
contact force by effecting a linear interpolation for a finite 
element modeling process. 

26. (NEW) A method according to claim 11, characterized in 
that barycentric coordinates are used to distribute the 
displacements and the forces of the points of application of the 
contact force between the extreme points of application of the 
contact force by effecting a linear interpolation for a finite 
element modeling process. 

27. (NEW) A method according to claim 9, characterized in that 
when the points of application of the contact forces between two 
objects in collision have been determined, during the step d) 

-15- 

WEINGARTEN, SCHURGIN, 
GAGNEBIN & LEBOVICI LLP 
TEL. <617) 542-2290 
FAX. (617) 451-0313 



Attorney Docket No. BDL-506XX 
Filed: Herewith 
TC Art Unit: 
Confirmation No. : 

the mechanical characteristics of the objects are transferred 
into the defined contact space in which the whole of a group of 
in contacts with n objects is processed, where it? and n are 
integers . 
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